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ABSTRACT 

Following the successful mapping of the last scattering surface by WMAP and 
balloon experiments, the epoch of the first stars, when Population III stars formed, is 
emerging as the next cosmological frontier. It is not clear what these stars' properties 
were, when they formed or how long their era lasted before leading to the stars 
and galaxies we see today. We show that these questions can be answered with 
the current and future measurements of the near-IR cosmic infrared background 
(CIB). Theoretical arguments suggest that Population III stars were very massive and 
short-lived stars that formed at z ~ 10 — 20 at rare peaks of the density field in the 
cold-dark-matter Universe. Because Population III stars probably formed individually 
in small mini-halos, they are not directly accessible to current telescopic studies. We 
show that these stars left a strong and measurable signature via their contribution 
to the CIB anisotropics for a wide range of their formation scenarios. The excess in 
the recently measured near-IR CIB anisotropies over that from normal galaxies can 
be explained by contribution from early Population III stars. These results imply 
that Population III were indeed very massive stars and their epoch started at ^; ~ 20 
and lasted past ztl3. We show the importance of accurately measuring the CIB 
anisotropies produced by Population III with future space-based missions. 

Subject headings: cosmology: theory - cosmology: observations - diff'use radiation - 
large scale structure of Universe 
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1. Introduction 

The cosmic infrared background (CIB) arises from the accumulated emission of galaxy 
populations spanning a large range of redshifts. The earliest epoch for the production of this 
background occurred when star formation first began, and contributions to the CIB continue 
through the present epoch. The CIB is thus an integrated summary of the collective star forming 
events, star-burst activity, and other luminous events in cosmic history to the present time. 

With the WMAP measurements (Bennett et al 2003) the structure of the last scattering 
surface has been mapped and a firm cosmological model has now emerged: the Universe is 
flat, dominated by the vacuum energy or an exotic quintessence field and is consistent with 
the inflationary paradigm and a cold-dark-matter (CDM) model. In this paper, we adopt 
the ACDM model with cosmological parameters from the WMAP and other observations: 
^^baryon = 0.044, h = 0.71, Q.rn = 0.3, Oa = 0.7, cTs = 0.84. Following the last scattering at z ^ 1000, 
the Universe entered an era known as the "Dark Ages", which ended with the star formation 
which produced Population III stars. The WMAP polarization results (Kogut et al 2003) show 
that the Universe had an optical depth since last scattering of r ~ 0.2, indicating an unexpectedly 
early epoch of the first star formation (z* ~ 20). Prom the opposite direction in z, optical and IR 
telescopes are now making progress into understanding the luminosity history during the most 
recent epoch of the universe (z < 5), but the period from recombination to the redshift of the 
galaxies in the Hubble Deep Field (HDP) remains largely an unexplored era. The first objects in 
the Universe probably formed in mini- halos of ~ 10^-^Mo (e.g. Miralda-Escude 2003) (barring the 
possible radiative feedback effects) , which would remain largely inaccessible to current instruments 
and conventional methods. 

A consensus based on recent numerical investigations is now emerging that fragmentation of 
the first collapsing clouds at redshift was very inefficient so that the first metal-free stars were 

extremely massive with mass IOOMq (Abel et al 2002, Bromm et al 1999). Such stars would 
live only a few million years, which is much less than the age of the Universe (~ 2 x 10^ years 
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at z = 20), making their direct detection still more difficult. It is unclear how long the era of 
Population III stars lasted and when the Population II stars formed with metallicities ~ 10~^ 
solar (Schneider et al 2002). On the other hand the net radiation produced by these massive stars 
will give substantial contributions to the total diffuse background light (Rees 1978) and since their 
light is red-shifted much of that contribution will be today in the infrared bands (Bond et al 1986) 
contributing significantly to the near-IR CIB. 

Observationally, the CIB is difficult to distinguish from the generally brighter foregrounds 
from the local matter within the Solar system, and the stars and ISM of the Galaxy. A number 
of investigations have attempted to extract the isotropic component (mean level) of the CIB 
from ground- and satellite-based data (see Hauser & Dwek 2001 for a recent review). These 
analyses of the COBE data have revealed the CIB at far-IR wavelengths A > 100/im(Hauser et al 
1998), and probably at near-IR wavelengths from 1-3 /xm (Dwek k, Arendt 1998, Gorjian et al 
2001, Cambresy et al 2001) with additional support from analysis of data from IRTS (Matsumoto 
et al 2000). However, none of the reported detections of the isotropic CIB are very robust 
(especially at near-IR wavelengths), because all are dominated by the systematic uncertainties 
associated with the modeling and removal of the strong foreground emission. Furthermore, the 
near-IR colors of the mean CIB do not differ greatly from those of the foregrounds, limiting the 
use of spectral information in distinguishing the true CIB from residual Galactic or solar system 
emission. Because of the difficulty of exactly accounting for the contributions of these bright 
foregrounds when measuring the CIB directly, and the difficulty in detecting all the contributing 
sources individually, it was proposed to measure the structure or anisotropy of the CIB via its 
angular power spectrum (Kashlinsky et al 1996). For a relatively conservative set of assumptions 
about clustering of distant galaxies, fluctuations in the brightness of the CIB have a distinct 
spectral and spatial signal, and these signals can be more readily discerned than the actual mean 
level of the CIB. 

In Fig. 1, we summarize the current near-IR CIB measurements and compare the measured 
CIB at 1.25 and 2.2 /xm to the total flux from the observed galaxy populations. There is a 
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Fig. 1. — Mean levels of the near-IR CIB from IRTS (crosses from Matsumoto et al 2000 and from COBE DIRBE (filled 
diamonds - at 2.2 nm from Gorjian et al 2001 and at 3.5 /im from Dwek & Arendt 1998. Cumulative flux from observed J and 
K band galaxies is shown with the filled circles (Cambresy et al 2001, Gardner 1996, Kashlinsky & Odenwald 2000b) Open 

triangle shows the CIB results of the DIRBE J-band analysis from Cambresy et al 2001. The insert shows the cumulative flux 
distribution in nW/m-^/sr as a function of K magnitude from galaxies from deep galaxy surveys (Totani et al 2001). 

statistically significant excess over the CIB expected from normal galaxies; e.g. at 1.25 ^m (J 
band) and 2.2 /im (K band), while an integration of the measured galaxy counts gives only 30-40% 
of the observed CIB (Kashlinsky & Odenwald 2000b, Gardner 1996, Madau & Pozzetti 2000). If 
this deficit CIB does not come from undetected low surface brightness emission associated with 
local galaxies, it must come from fainter galaxies at high z. 

The measurements of CIB structure or fluctuations lead to similar conclusions (Kashlinsky 
& Odenwald 2000a, Matsumoto et al 2000, Kashhnsky et al 2002, Odenwald et al 2003). The 
detection of CIB fluctuations at ^ 0.7° from COBE DIRBE bands at 1.25 to 5 micron (Kashlinsky 
&: Odenwald 2000a) shows a signal 2 — 3 times that expected from normal galaxy populations 
evolving via a present day IMF with the observed star formation rates out to z ~ 3 — 4 (Jimenez &; 
Kashlinsky 1999). Recently, using combined 2MASS calibration data with long net exposures, it 
became possible to remove resolved galaxies and measure the near-IR CIB fluctuations spectrum 
from a few arcsec to a few arcmin. We argued that this signal originated from faint galaxies 



K band 
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(K>19) at high z (Kashhnsky et al 2002, Odenwald et al 2003), and is significantly larger that 
what is expected from normal galaxy populations. This suggests that the excess originates from 
still earlier epochs. 

2. CIB flux and Population III 

The near-IR CIB signals detected in the 2MASS long integration data (Kashlinsky et 
al 2002, Odenwald et al 2003) could come from the Population III stars (Magliochetti et al 
2003, Kashlinsky et al 1999, Cooray, A. et al 2003, Santos et al 2002, Salvaterra & Ferrara 2003). 
The argument can be generalized and shown to be nearly model-independent provided Population 
III stars were very massive. 

Massive, metal-free stars will be dominated by radiation pressure and would radiate close 
to the Eddington limit: L ~ i^Edd = ^^!^^M ~ 1.3 x lO^^M/M© erg/sec, where or is the 
cross section due to electron (Thompson) scattering. The energy spectrum for emission from 
these stars will be nearly featureless and close to that of a black body at T ~ 10^-^""^K (Schaerer 
2002, Tumlinson et al 2003). Unlike stars with higher metallicity. Population III stars are not 
expected to have significant mass-loss during their lifetime (Baraffe et al 2001). If sufficiently 
massive (^240Mo, see Bond et al 1984) such stars would also avoid SN explosions and collapse 
directly to black holes. In this case their numbers could significantly exceed the minimum that 
would be required to produce the metallicities observed in Population II stars. The lifetime 
of these stars will be independent of mass: tL 0^ eMc^/L ~ 3 X 10^ years, where e = 0.007 is 
the efficiency of the hydrogen burning. These numbers are in good agreement with detailed 
computations (Schaerer 2002). 

For a flat Universe the co-moving volume per unit solid angle contained in the cosmic time 
interval dt is dV = c(l -|- z)~^d'j^dt, where di is the luminosity distance. The flux per unit 
frequency interval from each of the Population III stars will be Lbi,/{1 + 2;)/47rd|, where bi, is 
the fraction of the total energy spectrum emitted per unit frequency and = u{l + z) is the 



-7- 



rest-frame frequency. (The Population III SED is normalized so that / bj^dv = 1). The co-moving 
mass density in these stars is / Mn{M, t)dM = i^baryong;;^/*, where /* is the fraction of the total 
baryonic mass in the Universe locked in Population III stars at time t. The net flux per unit 
frequency from a population of such stars with mass function n{M)dM is given by: 

d fLn{M,t)dM ^ , n ^ dV 

J^h = (1 + Z) (6.,) — = ^(6.,)(-)/.pba.yonic (1) 



Here {hi,) = J Ln{M, t)b„dM/ J Ln{M, t)dM denotes the mean Population III SED averaged over 
their initial mass function and {^) = J Ln{M,t)dM/ J Mn{M,t)dM. For a Gaussian density 
field /* ~ 5 X 10~^ — 3 X 10^'^, if on average Population III formed in 2-3 sigma regions (see 
below). Provided that L/M = constant, this result does not depend on the details of the initial 
mass function of Population III stars (cf . Fig. 5 in ref. Salvaterra &; Ferrara 2003) . (Note that for 
the present day stellar populations their mass-to-light ratio depends strongly on stellar mass and 
{^) is much smaller than the Population III value of A-KGrnpc/ar — 3-3 x IO^Lq/Mq leading to 
substantially smaller net fluxes). Assuming no significant mass loss (Baraffe et al 2001) during 
their lifetime ti for Population III stars, these stars will produce GIB flux of amplitude: 

3 1 ^ 1 

Note that c^/G is the maximal luminosity that can be achieved by gravitational processes 
and enters here because the nuclear burning of stars evolves in gravitational equilibrium 
with the (radiation) pressure. This result, eq. 2, has a simple meaning illustrating its 
model-independence: the cumulative bolometric flux produced is Lmax = c^/G ~ 10^^ Lq 
distributed over the surface of the Hubble radius, 4TrR%, times the model-dependent dimensionless 
factors. The spectral distribution of the produced diff'usc background will be determined by 
(biy) = J /*(1 -I- z)~^{bi,)dtl J /*(1 -I- z)~^dt, which is the mean SED averaged over the flrst star 
epoch (hereafter assumed to have duration t*). The value of i^max/47ri?|f is ~ 3 x 10^ nW/m^/sr so 
even with small values of ^^baryom e, /*, the net flux from Population HI stars would be substantial. 
The K-band roughly coincides with the Lyman limit (912 A) at ~ 20 so much of the emitted 
flux would contribute today in the near IR bands. If there was a population of massive black 
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holes, as seems likely to result from Population III evolution, their accretion processes will emit 
radiation at a much higher efficiency of e ~ 0.1 (and perhaps as high as e ~ 0.4 very close to the 
event horizon). At the same time, the K-band CIB excess over the total contribution from the 
observed galaxy populations (~ 9nW/m^/sr up to K~ 24), is ~ 10 — 15nW/m^/sr and can be 
completely accounted for by the contribution from Population III stars at high redshifts. 

The first stars formed in rare regions of the density field that reached the turn-around 
over-density while the bulk of the matter was still in a linear regime (density contrast < 1). 
These regions stopped expanding and began collapsing, forming the first generation of stars if 
certain conditions are met. For a Gaussian density field, as required for the CDM models, this 
would specify the value of at each time. The power spectrum of the mass distribution on large 
scales preserved the Harrison-Zcldovich power spectrum index (n ~ 1), but on scales less than 
the horizon scale at the matter-radiation equality, the spectrum was modified by the diflFerent 
growth-rates of sub- and super-horizon modes. This led to a gradual decrease of n toward the 
very small scale value of — 2, suppressing the small scale power and leading to late collapse 
and galaxy formation. In this picture, it is surprising that the the first stars started forming at 

as indicated by the WMAP polarization measurements, but the most straightforward 
explanation is that the earliest stars formed out of very high peaks of the density field, leading to 
smaller 

Assuming a spherical model for growth and collapse of density fluctuations, any density 
fluctuation that reached the dimensionless density contrast 6coi = 1-68 will collapse. Consequently, 
if (Tm is the rms fluctuation at z over a sphere containing the total mass M, then 77 = Scoi/ctm 
measures the numbers of standard deviations for that mass to collapse at that redshift. Here, 
(Jm is related to the underlying power spectrum via = ^ J P{k; z)W{krM)k'^dk where 
rM = 1-4:{M /lO^'^ h^^ MqY'''^ {i}/0.3)^^^^ h^^Mpc is the co-moving linear scale containing mass 
M, and W{x) = {3ji{x) /x)"^ . As they collapse, the primordial clouds will quickly heat up to 
their virial temperature, Tyir, and only efficient cooling will enable isothermal collapse, where gas 
pressure is smaller than gravity (or T < Tyir). 
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The details of what determines the collapse and formation of Population III stars are not yet 
clear, but molecular cooling is critical in the initial stages of metal-free collapse. Molecular cooling 
is effective at Tt^OO K (Santos et al 2002), but some simulations suggest that rotational H2 
cooling can effectively dissipate binding energy of the cloud only if T > 2000 K (Miralda-Escude 
2003). We have evaluated the rms mass density fluctuation, gm, for the ACDM model as a 
function of its virial temperature, assuming a molecular weight of /x = 1. The left panel of Fig. 2 
shows the resulting rj ys z for Tvir = 400K and 2000K. For a Gaussian density field and 77 > 1 the 
value of f^, ~ erfc(r7/\/2). One can sec that Population III stars had to form out of 2 to 3 sigma 
fluctuations, which for a Gaussian distribution correspond to /* varying between ~ 5 x 10~^ and 
~ 3 X 10~^. Departures from spherical symmetry would accelerate the collapse, decreasing 77 and 
increasing f^. Because rj is a decreasing function of decreasing z, and /* is a (rapidly) increasing 
function of decreasing rj, the average /* will be dominated by the late times of Population III 
evolution. 




Fig. 2. — Left: The number of standard deviations that correspond to Population III regions collapsing at z for ACDM 
power spectrum of the primordial density field. Triangles correspond to Tvir = 400K, asterisks to Tvir = 2 x lO^K. Right: 
cumulative bolometric flux by Population III stars forming at 2» = 20, and lasting until Zf shown on the horizontal axis. Solid 
hne corresponds to Tvir = 400K and dotted Une to Tvir = 2 x lO^K. Thick dashed line shows the total CIB flux from DIRBE 
measurements at J,K and L bands; the shaded region shows the uncertainty. The values adopted were 22.9 ± 7.0 KJy/sr in J 
band (Cambresy et al 2001) and 16.1 ± 4.4, 13.4 ± 3.7 KJy/sr in K, L bands (Wright 2003). 

The right panel of Fig. 2 shows that the total diffuse background light produced during the 
first star era should be substantial and measurable. Much of the diffuse flux produced by the 
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Population III stars will today be observable in the near-IR, and could be responsible for the 
observed excess near-IR CIB. If Zf is large the net bolometric flux from Population III is smaller, 
but because of the redshift effects a larger fraction of the bolometric flux will be observed today in 
the near-IR bands. Furthermore, the near-IR part of the CIB emission from Population III would 
be increased because much of the shorter wavelengths emission gets absorbed and re-emitted at 
these bands by Lyman and free-free emission from the ionized nebulae around Population III 
stars (Santos et al 2002). Dashed line shows the total CIB flux evaluated from the DIRBE 1.25 
to 3.5 //m bands. Taken at face value that would imply that the first star era lasted until Zf^ 15. 
However, Population III stars only contribute a fraction of the observed CIB. 



The contribution of Population III stars to the CIB will possess a spatial structure that is 
related to that of density fluctuations from which the stars formed. At particular spatial scales and 
wavelengths it may prove easier to distinguish the Population III contribution from the various 
foregrounds through the spatial structure rather than the measurement of the mean (isotropic) 
intensity. 

In the limit of small angles the 2-dimensional angular power spectrum, P2{q), of the CIB 
fluctuations that arise from clustering of Population III systems with the 3-dimensional power 
spectrum P3{k) is given by the modified Limber equation (Kashlinsky & Odenwald 2000a): 



where c/a is the comoving angular diameter distance. The power spectrum P3 is in turn related 
to the underlying ACDM model power spectrum and its evolution and growth. On larger angular 
scales {0t a few arcmin) the power spectrum is in the linear regime and we evaluate the growth 
rate of fluctuations using linear theory for the A-dominated cosmology. 

On sufficiently small scales the primordial power spectrum of the density field will be modified 



3. CIB anisotropies from Population III 




(3) 
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by non-linear gravitational effects. Fig. 3 shows the density field for the ACDM model at various 
high redshifts. The linear power spectrum was taken to be that of the ACDM for the WMAP 
cosmological parameters. The non-linear evolution of the density field was modeled analytically 
using the approximation from ref. (Peacock & Dodds 1996). For reference 1 arcmin subtends the 
comoving scale of ~ 1.5/i^^Mpc at z between 6 and 15 for WMAP cosmological parameters and 
a dark energy equation of state with w = —\. The dotted line shows the linearly evolved density 
field at z = 15. On sub-arcminute scales the density field is in the quasi-linear to non-linear regime 
(density contrast iiO.2) for the ACDM model. There the spectrum due to clustering evolution 
was modified significantly and the fluctuations amplitude has increased from its primordial value, 
especially since the effective spectral index on these scales for ACDM model was — 2. 

100.000 p N 




0.001 r 1 

0.001 0.010 0.100 1.000 10.000 100.000 
k"' (h^'Mpc) 

Fig. 3. ACDM density field at z = 6, 10, 15 (red, green and blue). Dotted line corresponds to linearly evolving density 

field at z = 15. 

To within a factor of order unity, the square of the fractional fluctuation of the CIB on angular 
scale ~ ir/q is Sqj^ = {{Sly)"^) / ~ IjJ^q'^P2{q)/2iT. The meaning of eq. 3 can be demonstrated 
by assuming (i/,^//(it=constant during the Population III phase of duration t^,. In this case the 
fractional fluctuation due to clustering of early star systems becomes: 

5l,^ = ^J /\\qd-/-z)dt (4) 

Here A(A;) = ^ J^3(fc) ji/2 -g ^j^g fluctuation in number of sources within a volume k~'^ct^. In 
other words, the fractional fluctuation on angular scale vr/g in the CIB from Population III stars 
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is given by the average value of the r.m.s. fluctuation from Population III spatial clustering over 
a cylinder of length ct* and diameter ~ At z = 20 the Universe is ~ 2 x 10^ years old 

which is much larger than the age of the individual Population III stars, t/, ~ 3 x 10^ years. If 
the first stars lasted for only one generation forming at z*, the relative CIB fluctuations will be 
SciB ~ ^/^^A{qd'2^■,z)\z,. 

Several points are worth noting about eq. 3: Firstly, the value of A{k) is inversely proportional 
to -y/U and thus measures the duration of the first star era. The density perturbations grow 
with decreasing z, so most of the contribution to the integral in eq. 3 comes from the low end of 
z. The overall dependence on wins out and the longer the Population III phase, the smaller 
are the relative fluctuations of the CIB from them. Secondly, in the Harrison-Zeldovich regime 
of the power spectrum, P3 oc k, one would have Sqib q^'^. Finally, the transition to the 
Harrison-Zeldovich regime occurs in the linear regime at all relevant redshifts and happens at the 
co-moving scale equal to the horizon scale at the matter-radiation equality. This feature would 
occur at different angular scales as the redshift of the Population III era decreases. All this would 
allow probes of the epoch of the first stars, its duration, and the primordial power spectrum at 
high redshifts on scales that are currently not probed well. Interestingly, a short duration for the 
first stars will lead to smaller CIB flux, but larger relative fluctuations and vice versa. 

In addition to the small angular scale increase due to non-linear gravitational evolution, the 
fluctuations in the distribution of Population III systems will be amplifled because, as Fig. 2 
shows, these systems had to form out of rare peaks of the primordial density field. Consequently, 
their 2-point correlation function will be amplifled (Kaiser 1984, Jensen &; Szalay 1986, Politzer 
& Wise 1984, Kashlinsky 1987, Kashlinsky 1998) over that given by the ACDM power spectrum, 
= ^ J Ps{k)jo{kr)k'^dk. Adopting a Press-Schechtcr prescription for the formation of collapsed 
objects (Press & Schechter 1974) from a correlated density field leads to the following expression 
for the amplification of the 2-point correlation function (Kashlinsky 1998): 
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where are Hermite polynomials. At large values of r) this reduces to ~ exp{r)'^^/S'^^i) — 1. 
At small values of ^ <C 1 (and sufficiently large angular scales) the amplification factor is linear 
and equal to ~ 'h'^/^coV very small scales the amplification is non- linear in ^ and will exceed 
the linear amplification value leading to larger CIB fluctuations. We computed the resultant CIB 
fluctuations using the linear bias amplification approximation, but note that this gives a lower 
limit on the resultant CIB fluctuations. 

1000 F — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 1 




5 10 15 20 

Zf 

Fig. 4. — Mean projected surface density of Population III systems along the line of sight, plotted vs the redshift marking 
the end of the Population III era. Solid and dotted lines correspond to primordial clouds containing Population III with Tvir > 
400 K and 2000 K respectively. 

There will also be shot-noise fluctuations due to individual Population III systems entering 

— 1/2 

the beam. The relative magnitude of these fluctuations will be A^b^ajn' where A^bcam is the number 
of the Population III systems within the beam. This component will be important at very small 
angular scales, where Ab 

earn 1) £md will Contribute to the power spectrum: .PsN — l/'^2) where 
7X2 = c / JT-ac^iCl + z)''^dt is the projected angular number density of Population III systems. 
Detection of the shot-noise component in the CIB power spectrum at small angular scales will give 
a direct measure of both the duration of the first star era and constrain the makeup and masses 
of the Population III systems. The three-dimensional number density can be evaluated from the 
Press-Schechter (Press & Schechter 1974) formalism, but in the appropriate limits [r] > 1) it 
becomes ns ~ /9baryon-^~''^erfc(?7(z)/\/2). Fig. 4 shows that the mean projected surface density of 
Population III systems n2 for Tyir > 400 K and 2000 K. One can see that unless the Population 
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III era was very short, the shot noise correction to the CIB would be smaU on scales greater than 
a few arcseconds. 

Essentially all of the flux emitted at rest wavelength less than 912 A will be absorbed by 
the Lyman continuum absorption in the surrounding medium (Yoshii &; Peterson 1994) and 
the contribution to the CIB will be cut off at A < 912(1 + Zf) A. At longer wavelengths the 
flux, including emission from the nebulae around each star, will propagate without significant 
attenuation. Population III stars emit as black bodies at log T ~ 4.8 — 5 (Schaerer 2002) and the 
Lyman limit at ^; ~ 20 is shifted to ~ 2/im in the observer's frame. The 2MASS J band filter 
contains emission out to 1.4 fim and represents the shortest wavelength where excess in the CIB 
over that from normal galaxies has been measured (Cambresy et al 2001, Kashlinsky et al 2002). If 
the measured excesses in CIB fluctuations (Kashlinsky et al 2002, Kashlinsky &: Odcnwald 2000a) 
and isotropic component (Cambresy et al 2001) at J band are indeed attributable to Population 
III, their era must have lasted until 14. 

Fig. 5 shows the resultant CIB fluctuations from Population III stars at 1.25, 2.2 and 3.5 
micron or J, K and L photometric bands. The assumed CIB fluxes from Population III are listed 
near the tops of the panels. At 2.2 micron the galaxies contribute around 9 nW/m^/sr whereas 
the net CIB flux is around 20-25 nW/m^/sr ; hence the assumed excess of 10-15 nW/m^/sr is 
attributed to the Population III stars for the piuposes of this discussion and is taken to be 12 
nW/m^/sr in the figure. The situation is similar at 3.5 micron (Dwek Sz Arendt 1998, Gorjian et 
al 2001) and the excess was taken to be 5 nW/m^/sr (Wright 2003). The excess in the J band 
over the total flux from normal galaxies was taken from (Cambresy et al 2001). The clustered 
component of the CIB fluctuation displayed in the figure can be rescaled to other CIB fluxes via 
eq. 4. Thin colored lines assume constant rate of production of the Population III fraction of the 
CIB and the flat SED with = constant. All Population III stars were assumed to start forming 
at = 20, but different colored lines correspond to different values oi Zf, the redshift of the end 
of the first star era. Thick colored lines correspond to a Population III contribution produced by 
black-body emitters in the Rayleigh- Jeans regime {bi, oc u^) with the rate of flux production oc /*. 
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We also computed CIB fluctuations with the rate of flux production oc for more realistic SED 
from Santos et al (2002) which include the Lyman-a and free-free emission from the surrounding 
nebula. Thin and thick dashed colored lines correspond to two possible extremes: SEDs from 
their Fig. 3 (completely opaque nebula to ionizing photons) and Fig. 5 (completely transparent 
nebula). Because they do not differ significantly from the more simple SEDs the lines are shown 
only for some values of to avoid overcrowding in the Figure. The graphs show the relative 
model-independence of the relative CIB fluctuations on the SED for a fixed total CIB flux; the 
latter contains the bulk of information on the details of Population III emission distribution 
between various NIR bands. For a given Population III CIB isotropic component, the fluctuations 
are largest for large values of Zj because the duration of the first star era is the smallest. Because 
of the Lyman continuum absorption there would be no appreciable CIB fluctuations at J band if 
z/i; 14 (no blue lines in the left panel) and the contribution to the other scenarios will come from 
z ~ 13, rather than z = 20 when the first star era was assumed to start. This would be broadly 
consistent with the limits on sky brightness fluctuations at 4 fim from Xu et al (2003). 

Both DIRBE and 2MASS data indicate CIB anisotropics at amplitudes larger than the 
contribution from normal galaxy populations (Kashlinsky & Odenwald 2000a, Kashlinsky et al 
2002) and are consistent with significant contributions due to Population III. However, because 
the measured signal contains contributions from remaining galaxies (all galaxies for DIRBE and 
K > 19 galaxies for 2MASS), it is difficult to isolate the contribution from the Population III 
stars. 

4. Prospects for measuring CIB from Population III 

Population III stars, if massive, should have left a unique and measurable signature in the 
near-IR CIB anisotropics over angular scales from ^ 1 arcminute to several degrees as Fig. 5 
shows. As is the case for measurement of the mean isotropic CIB, detection of these fluctuations 
depends on the identification and removal of various foreground emission (and noise) contributions: 
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atmosphere (for ground-based measurements), zodiacal light from the Solar system, Galactic 
cirrus emission, and instrument noise. 

The dashed line in Fig. 5 shows the small scale atmospheric fluctuations at 2 after 
one hour of integration of ground based measurement (Odenwald et al 2003). At large angular 
scales atmospheric gradients will become important, making measurements there even more 
difficult. Atmospheric effects can be highly variable on a wide range of time scales, yet they 
can be completely avoided with space-based experiments. Zodiacal emission from interplanetary 
dust (IPD) is the brightest foreground at most IR wavelengths over most of the sky. There 
are some structures in this emission associated with particular asteroid families, comets, and 
an earth-resonant ring, but these structures tend to be confined to low ecliptic latitudes or 
otherwise localized. The main IPD cloud is generally modeled with a smooth density distribution. 
Observationally, intensity fluctuations of the main IPD cloud have been limited to < 0.2% at 
25 /um (Abraham et al 1997). Extrapolating this limit to other wavelengths using the mean 
high-latitude zodiacal light spectrum yields the limits shown in Figure 5. Because the Earth is 
moving with respect to (i.e. orbiting within) the IPD cloud, the zodiacal light varies over time. 
Likewise, any zodiacal light fluctuations will not remain fixed in celestial coordinates. Therefore 
repeated observations of a field on timescales of weeks to months should be able to distinguish and 
reject any zodiacal light fluctuations from the invariant Galactic and CIB fluctuations. 

Intensity fluctuations of the Galactic foregrounds arc perhaps the most difficult to distinguish 
from those of the CIB. Stellar emission may exhibit structure from binaries, clusters and 
associations, and from large scale tidal streams ripped from past and present dwarf galaxy 
satellites of the Milky Way. At mid- to far-IR wavelengths, stellar emission is minimized by virtue 
of being far out on the Rayleigh-Jeans tail of the stellar spectrum (apart from certain rare classes 
of dusty stars). At near-IR wavelengths stellar emission is important, but with sufficient sensitivity 
and angular resolution most Galactic stellar emission, and related structure, can be resolved and 
removed. IR emission from the ISM (cirrus) is intrinsically diffuse and cannot be resolved. Cirrus 
emission is known to extend to wavelengths as short as 3 fim. Statistically, the structure of the 
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cirrus emission can be modeled with power-law distributions. Using the mean cirrus spectrum, 
measurements made in the far-IR can be scaled to 3.5 |Ltm, providing the estimated fluctuation 
contribution from cirrus that is shown in Figure 5. The extrapolation to shorter wavelengths is 
highly uncertain, because cirrus (diffuse ISM) emission has not been detected at these wavelengths, 
and the effects of extinction may become more significant than those of emission, but the cirrus 
contribution is expected to be several orders of magnitude lower than the CIB fluctuations. 

The current detections of CIB anisotropics from DIRBE and 2MASS are above all these 
foregrounds. The level of anisotropy compares favorably with theoretical expectations and may 
indeed contain a signiflcant signal from Population III stars. However, the anisotropy may also 
contain contributions from galaxies made up of Populations II and I stars, which may be signiflcant 
and uncertain at high z because of unknown evolution effects. Hence it is important to be able 
to remove galaxy contribution in future measurements. Galaxies with normal stellar populations 
fainter than K ~ 24 contribute little to the CIB. At fainter magnitudes (and larger redshifts) 
the emission is likely to be dominated by Population HI stellar populations and protogalaxies. 
Because each unit mass of Population HI emits a factor of ~ 10^ more light that the normal 
(Population I, II) stellar populations, the diffuse flux from (massive) Population III stars would 
dominate emission from faint normal galaxies (say K > 24) . The figure shows that CIB fiuctuations 
from Population HI would be the dominant source of diffuse light fluctuations on arcminute and 
degree scales even if Population HI stars epoch was briefer, and their diffuse flux smaller, than 
the current CIB numbers suggest. Their angular power spectrum should be very different from 
other sources of diffuse emission and its measurement thus presents a way to actually discover 
Population HI and measure the duration of their era and their spatial distribution. The latter 
would provide direct information on primordial power spectrum on scales and at epochs that are 
not easily attainable by conventional surveys. This measurement would be imperative to make 
and is feasible with the present day space technology. 
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5. Observational requirements 



There are several ways in which observations could verify the analysis of this paper. The 
James Webb Space Telescope (JWST) will be capable of resolving individual Population III objects 
in their supernova phase, or as young black holes accreting substantial mass flows. Alternatively, 
if the Population III objects are sufficiently clustered in space and time, the aggregations would 
resemble proto-galaxies and could be recognized. Observations of these objects are among the 
foremost objectives of the JWST. 

The following would be the optimal set of parameters for an observation designed to map the 
contribution of the first star era to GIB anisotropies: 1) In order to make certain that the signal is 
not contaminated by distant galaxies with normal stellar populations, one would need to conduct 
a deep enough survey in order to identify and eliminate normal galaxies from the field. In practice, 
as Fig. 1 indicates, going to K~ 24 would be sufficient. With the current technology this would 
take 1.5 hours for a 1 meter space telescope. 2) A direct signature of Population III signal is that 
there should be no GIB fiuctuations at wavelengths < 912(1 + Zf)~^ A. That spectral drop would 
provide an indication of the epoch corresponding to the end of first star era. At longer wavelengths 
the spectral energy distribution of the GIB fiuctuations would probe the history of energy emission 
and its re-distribution by the nearby gas during the first star era. At wavelengths 1 10/^m zodiacal 
light fluctuations may become dominant. 3) On angular scales from a few arcminutes to ~ 10° 
Population III would produce GIB anisotropies with a distinct and measurable angular spectrum, 
but its measurement will be limited by the cosmic or sampling variance. If the power spectrum 
is determined from fraction /sky of the sky by sampling in concentric rings of radius q ^ tt/O and 
width Aq, the relative uncertainty on P2{q) will be ~ WW fskJ"^- order to get 



accurate and independent measurements, it is necessary to have small Aq/q ~ 0.05 — 0.1; thus for 
reliable results on scales up to 9 one would need to cover area a few times larger. Shaded regions 
in the left panel of Fig. 5 show the one sigma cosmic variance uncertainty with Aq/q = 0.05 for a 
total of 1, 100 and 1,000 deg^ areas covered. In order to get reliable results on scales up to 9^1° 
one would need to observe areas of 10° across. 
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The black box in Fig. 6 shows what parameters the ideal space mission should have for such 
a measurement with cross-hatched region in the right panel displaying the magnitude range where 
Population III stars are expected to dominate CIB anisotropics . For comparison other currently 
operating or planned NASA missions in the near-IR are also shown. Some of the datasets overlap 
with the required parameter space, but a large part of the desired parameter space can not be 
reached with these other data sets. 

Direct observations of the CIB fluctuations with the needed fidelity and angular range is not 
easy. The foreground objects (stars and low-redshift galaxies, and the interstellar medium) are 
bright and abundant. However, they do not have the same spatial and spectral distributions as 
the predicted CIB fluctuations. The first observational requirement is to measure the fluctuations 
as a function of spatial scale and wavelength, with sufficient corroborating evidence to show that 
the fluctuations are not the residual of some foreground process. 

The predicted fluctuations turn down at angular scales greater than subtended by the horizon 
at decoupling seen at the first star era, and this must be verified. Masking out contributions 
from individual galaxies requires angular resolutions of the order of 1 arcsec and adequate depth 
to remove galaxies of K < 23 — 24. Observing fluctuations on the order of 3 degrees in size with 
reasonable resolution of the spatial scale and reasonable signal-to-noise ratio requires a map of the 
order of 1000 deg^ in size, about 1/40 of the whole sky. Taken at face value, we need a map of 
the order of ~ 10^° pixels, to be built up as a mosaic of exposures with individual detectors and 
telescope pointings. 

There are many technical challenges associated with such a project. Ideally the observations 
would be taken in space, to avoid the limitations of atmospheric noise, and to enable observations 
at the full range of wavelengths from 1 to 5 /xm. In addition, we would require software to 
stitch together exposures from different detectors, taken at different times with different telescope 
orientations, and different amounts of zodiacal light and cosmic ray hits on the detectors, while 
preserving accurate calibration. It is essential to preserve large-scale photometric accuracy to 
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enable measurements of fluctuations on scales larger than a single field of view (Fixsen et al 2000) . 
To understand and eliminate the effects of interstellar dust and PAH clouds, it will be necessary 
to provide photometric bands that maximize the detectability of these sources of interference. 
In addition, observations over wide ranges of foreground brightness would enable recognition 
of unexpected effects. Such a project would produce a large database of space observations of 
low-contrast, extended features, and could lead to unexpected discoveries. 

To summarize, the measurement of CIB fluctuations from the Population III star era is 
imperative to make and is feasible with the present day space technology. 
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Fig. 5. CIB fluctuations, -y/ q'^P2/'2-K, from Population III stars forming at z» = 20. Continuous blue, green and red lines 

correspond to Zf = 15, 10,6 respectively. Thin colored solid lines correspond to the case of constant CIB flux rate production 
and by = const. Thick colored solid lines correspond to dl^/dt oc ftb^i with bu oc u^. Lower lines correspond to Tvir = 400K, 
and top lines to T^ir = 2000K. Colored dashed lines show the CIB fluctuations for more realistic SED from Santos et al (2002) 
which include the Lyman-Q and free-free emission from the surrounding nebula. They mostly overlap with the other SED 
models showing the relative model-independence of the results. Thin and thick dashed colored lines correspond to SED from 
Fig. 3 (completely opaque nebula to ionizing photons) and Fig. 5 (completely transparent nebula). Thin dashed lines show 
atmospheric fluctuations from ground based 2MASS measurements after 1 hour of exposure. Thick long dashes denote the upper 
limit on zodiacal light fluctuations from Abraham et al 1997 scaled to the corresponding band. Thick dashed-triple-dotted line 
denotes cirrus fluctuations: these are upper limits at 2.2 micron and an estimate from Kiss et al 2003 at 3.5 micron. In the 
middle panel dot-dashed lines correspond to shot noise from galaxies fainter than K=24 (thickest), K=23 and K=20 magnitudes 
(thin). The K-band faint galaxy counts data used were taken from Totani et al 2001. Diamonds with error bars show the CIB 
fluctuations at ^ 0.7° from the COBE DIRBE fluctuations analysis (Kashlinsky & Odenwald 2000a). Note that because of 
the large DIRBE beam, the results include contribution from all galaxies as well as other sources such as Population III. Dots 
with error bars show the K band CIB fluctuation from deep integration 2MASS data (Kashlinsky et al 2002. The 2MASS data 
shown in the figure were taken for the the patches for which galaxies were removed brighter than K ~ 19'" (Odenwald et al 
2003). The cosmic variance 1-sigma uncertainty is shown with shaded regions in the left panel. The darkest shade corresponds 
to a total of 1,000 deg^ observed, the intermediate shade corresponds to a total of 100 deg^ and the lightest shade corresponds 
to 1 deg-^. 
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Fig. 6. Black rectangular box shows the regions of the parameters for a mission to measure the CIB anisotropics 

from Population III thereby detecting the first star era. Cross-hatched region in the right panel displays the magnitude 
range where Population III stars are expected to dominate CIB anisotropies. Other colors correspond to the currently 
existing datasets or planned NASA missions useful for this purpose: red is the Hubble Deep Field from HST 
(http://www.stsci.edu/ftp/science/hdf/hdf.html), dark blue is for the deep 2MASS data (Kashlinsky et al 2002, Odenwald 
et al 2003), green is the NASA SIRTF GOODS project (http://www.stsci.edu/ftp/science/goods/), light blue corresponds to 
the planned NASA WISE MIDEX mission (http://www. astro. ucla.edu~wright/WISE), and orange color corresponds to the 
field-of-view of the NASA James Webb Space Telescope planned for the end of the decade (http://jwst.gsfc.nasa.gov). 
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